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^ <—j Abstract 

> \ 

In Internet environment, traffic flow to a link is typically modeled by superposition of 
ON/OFF based sources. During each ON-period for a particular source, packets arrive 
according to a Poisson process and packet sizes (hence service times) can be generally 
distributed. In this paper, we establish heavy traffic limit theorems to provide suitable 
approximations for the system under first-in first-out (FIFO) and work-conserving ser- 
. vice discipline, which state that, when the lengths of both ON- and OFF-periods are 

(T) | lightly tailed, the sequences of the scaled queue length and workload processes converge 

weakly to short-range dependent reflecting Gaussian processes, and when the lengths of 
. ON- and/or OFF-periods are heavily tailed with infinite variance, the sequences converge 

weakly to either reflecting fractional Brownian motions (FBMs) or certain type of long- 
range dependent reflecting Gaussian processes depending on the choice of scaling as the 
number of superposed sources tends to infinity. Moreover, the sequences exhibit a state 
space collapse- like property when the number of sources is large enough, which is a kind 
■ of extension of the well-known Little's law for M/M/l queueing system. Theory to jus- 

, tify the approximations is based on appropriate heavy traffic conditions which essentially 

mean that the service rate closely approaches the arrival rate when the number of input 
sources tends to infinity. 

Keywords: reflecting fractional Brownian motion, reflecting Gaussian process, long- 
range dependence, queueing process, weak convergence 
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1 Introduction 

ON/OFF sources are widely used to model voice, video and data traffics in telecommunication 
systems (see, e.g., Jain and Routhier [16], Nikolaidis and Akyildiz [23], Taqqu et al. [28], 
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Leland et al. [19], Paxson and Floyd [26]). In particular, stochastic modeling of queueing 
systems with ON/OFF long-range dependent data has become an active area of research. In 
contrast to most of the existing achievements in this field, which are based on fluid models 
whose outputs are deterministic with constant (service) rates and whose inputs are certain 
types of long-range dependent fluid sources (e.g. Norros [24], Massoulie and Simonian [21], 
Debicki and Mandjes [10], Debicki and Palmowski [11]), we will model our queueing system 
with general service time distribution and the input as a superposition of Poisson ON/OFF 
point processes to better capture the variation of packet sizes and the behavior of real packet 
traffic. Concretely, for a particular source, packets arrive according to a Poisson process 
during each ON-period. For such a source, the corresponding traffic exhibits long-range 
dependence (see, for instance, Ryu and Lowen [27]) when the lengths of ON- and/or OFF- 
periods are heavily tailed with infinite variance. Besides the assumption on the service time 
distribution, our system is further supposed to operate under FIFO and work-conserving 
discipline. 

A special case of the above queueing model is discussed in Cao and Ramanan [3], where 
the distributions of ON- and OFF-periods are assumed to be Pareto and exponential re- 
spectively and packet sizes are supposed to be constant. They show that the sequence of 
probabilities that steady state unfinished works exceed a threshold tend to the correspond- 
ing probability assuming Poisson input process when the number of input sources tends to 
infinity. Currently, it is not clear whether their result can be extended to the more general 
model as presented above. Furthermore, the dependence of the convergence rates on various 
parameters of the system is not shown in their result, e.g., the relationship between p N (traffic 
intensity, utilization level) and iV (the number of sources). 

Due to the above reasons, we will study our queueing system by employing some other 
method. Under heavy traffic conditions (suitable relationships between p and N such that 
the service rate closely approaches the arrival rate when N tends to infinity), we will show 
that, when the lengths of both ON- and OFF-periods are lightly tailed, the sequences of 
the scaled queue length and workload processes converge weakly to short-range dependent 
reflecting Gaussian processes, and when the lengths of ON- and/or OFF-periods are heav- 
ily tailed with infinite variance, the sequences converge weakly to either reflecting fractional 
Brownian motions (FBMs) or certain type of long-range dependent reflecting Gaussian pro- 
cesses depending on the choice of scaling as the number of input sources tends to infinity. 
Moreover, the sequences exhibit a state space collapse- like property when N is large enough, 
which is a kind of extension of the well-known Little's law for M/M/l queueing system. 

Our heavy traffic limits set up certain connection between the above physical queueing 
systems and some existing fluid queueing models. For example, for a fluid model with constant 
output rate and FBM input, the stationary queue content distribution is asymptotically 
Weibullian (e.g., Norros [24], Massoulie and Simonian [21], Husler and Piterbarg [15], and 
more generally, as summarized in Whitt [29]), namely, the probability of exceeding buffer 
level b is roughly of the form exp(— 6 2 ( 1- ^)) if FBM is characterized by Hurst parameter H. 
The result can be applied to derive corresponding probability for our reflecting FBM after 
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properly managing parameters. 

Concerning heavy traffic limit theorems for queueing systems with long-range dependent 
inputs, there are only a few achievements until now besides the one mentioned above in [3]. In 
Debicki and Mandjes [10] and Debicki and Palmowski [11], authors studied a fluid queueing 
system with constant output rate and a superposition of ON/OFF fluid input sources. In 
Konstantopoulos and Lin [18] and Majewski [20], instead of discussing superposition prob- 
lem, authors considered a single class and feedforward multiclass queueing networks with 
long-range dependent interarrival and service time sequences respectively. The current limit 
theorems are the supplements of these existing results. In justifying our reflecting FBM 
approximation, we will adopt the simultaneous limit regime related to FBM in Mikosch et 
al. [22], in which both N (the number of sources) and T (the time-scaling parameter) go to 
infinity at the same time. This procedure provides us some convenience in employing some 
ingredient developed in [22] to establish the weak convergence for our scaled queue length 
and workload processes. 

One last point we wish to mention is that we have employed our theorem on reflecting 
Gaussian processes in the current paper to provide a reasonable interpretation (in Dai [9]) to 
some well-known large-scale computer and statistical experiments conducted by Cao et al. [2], 
Cao and Ramanan [3], which reveal some gapes between their simulation findings and the 
existing theory on heavy-tail and long range dependence. In Dai [9], the author finds out that 
all the 'heavy-tail' random variables used in computer and network simulations are truncated 
versions of their real heavy-tail counterparts due to the limitations of computer hardware 
and softwares, and hence they are not heavily tailed ones. So, by combining the findings in 
Dai [9] and the theorem in the current paper, we claim in Dai [9] that the findings in Cao 
et al. [2], Cao and Ramanan [3] are more close to practice and but not to the mathematical 
assumptions imposed in their models since their simulations are computer-based ones. 

To be convenient for readers, here we summarize some frequently used notations and 
terminologies throughout the paper. First, we recall the definition of u.o.c. convergence. 
For a function / : [0, oo) — > R and t > 0, put = sup 0<s<t |/(s)|, then a sequence of 
functions f n : [0,oo) — > R is said to converge uniformly on compact sets (u.o.c.) to / if 
for each t > 0, ||/ ra — f\\t — > as n — > oo. Second, we use Cf,(R) to denote the set of all 
bounded and continuous functions / and C{R) to denote the set of all continuous functions 
over the real number space R, which are endowed with the uniform topology. Third, we use 
De[0, oo) to denote the Skorohod topological space, i.e., the space of E- valued functions that 
are right continuous and have left-hand limits, which is endowed with the Skorohod topology 
(see, e.g., Ethier an Kurtz [12], Billingsley [1]). Fourth, we use i.i.d to denote independent 
and identically distributed, use a.s. to denote almost surely, use '=>' to denote 'converge in 
distribution' or equivalently 'converge weakly', and use '~' to denote 'equals approximately'. 

The rest of this paper is organized as follows. In Section 2, we formulate our model, and 
in Section 3, we present our main theorems and they are proved in Section 4. 
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2 Queueing model formulation 



In this section, we consider a queueing system with general service time distribution and 
with N i.i.d. Poisson ON/OFF input sources. Concretely, a Poisson ON/OFF source n £ 
{1,..,N} consists of independent strictly alternating ON- and OFF-periods, moreover, it 
transmits packets to a server according to a Poisson process with interarrival time sequence 
{u n (i),i > 1} and rate A if it is ON and remains silent if it is OFF. The lengths of the ON- 
periods are identically distributed and so are the lengths of OFF-periods, and furthermore, 
both of their distributions can be heavily tailed with infinite variance. Specifically, for any 
distribution F, we denote by F = 1 — F the complementary (or right tail) distribution, and 
by F\ and F<i the distributions for ON- and OFF-periods with probability density functions 
fi and f2 respectively. Their means and variances are denoted by /U, and of for i = 1,2. In 
what follows, we assume that as x — > oo, 

(2.1) either Fi{x) ~ x~ ai Li(x) with 1 < a« < 2 or of < oo, 
where Lj > is a slowly varying function at infinity, that is, 

lim - - X } = 1 for any t > 0. 

Note that the mean Uj is always finite but the variance of is infinite when «j < 2, and 
furthermore, one distribution may have finite variance and the other has an infinite variance 
since F\ and F% are allowed to be different. The sizes of transmitted packets (service times) 
form an i.i.d. random sequence {v N (i) = v(i) /fi N ,i > 1}, where /j, n is the rate of transmission 
corresponding to each N and {v (i) : i > 1} is an i.i.d. random sequence with mean 1 and 
variance of, moreover, {v (i) : i > 1} is independent of the arrival processes. 

To derive our queueing dynamical equation, we introduce more notations. For a single 
source n e {1, N}, it follows from the explanation in Mikosch et al. [22] that the alternating 
ON/OFF periods can be described by a stationary binary process W n = {W n (t),t > 0}: 
W n (t) = 1 means that input traffic is in an ON-period at time t and W n {t) = means that 
input traffic is in an OFF-period, and moreover, the mean of W n is given by 

(2.2) 7 = EW n (t) = P(W n (t) = 1) = /^/(/i! + /i 2 ). 

Let T n (t) denote the cumulative amount of time which the nth source is ON during time 
interval [0,t], that is, 

(2.3) T n (t) = fw n (s)ds. 

Jo 

Let A n (t) be the total number of packets arrived at the server from the nth source during 
[0, t], namely, 

(2.4) A n {t) = sup jm,f> n (z) < T n (t) j , 
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which exhibits long range dependence if o~\ and 02 are not finite simultaneously (see, for in- 
stance, Ryu and Lowen [27]). Moreover, let A N (t) be the total number of packets transmitted 
to the server by time t summed over all N sources, that is, 

N 

(2.5) A N (t) = J2Mt), 

n=l 

and let S N (t) be the total number of packets that finished service at the server if her keep 
busy in [0, t], that is, 

(2.6) S N (t) = sup {m,V N (m) <t}, 
where 

m 

(2.7) V N (m) = Y,v N d)- 

Then the queue length process Q N (t) which is the number of packets including the one being 
served at the server at time t can be represented by 

(2.8) Q N {t) = A N {t) - S N (B N (t)), 

where we assume that the initial queue length is zero for convenience, B N (t) is the cumulative 
amount of time that the server is busy by time t. In the following analysis, we will employ 
FIFO and non-idling service discipline under which the server is never idle when there are 
packets waiting to be served. Hence the total busy time can be represented as 

B N (t)= f I{Q N (s) >0}ds, 
Jo 

where /{•} is the indicator function. Finally, we introduce the below workload process which 
measures the delay of a packet staying in the system, 

(2.9) L N (t) = V N {A N (t)) - B N (t). 



3 Heavy traffic limit theorems 

We are interested in the behaviors of the queueing process Q N (-) and the workload process 
L N (■) under suitable scaling and under the condition that the load of the server closely 
approaches the service capacity when the source number N gets large enough. In order 
to state our main theorems, we introduce the below notations for convenience, which are 
adapted from Taqqu et al. [28]. When 1 < ai < 2, set = (r(2 — aj))/(aj — 1). When 
of < oo, set cti = 2, Li = 1 and dj = of /2. Moreover, let 

b= lim 

x^oo Loix) 
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If < b < oo (implying a\ = a 2 and b = Hindoo L\{x) / L 2 {x)) , set a m i n = ot\, 

f oi\ 2 2{n 2 2 aib + n\a 2 ) 

(3.1) 7r = rzrrn r and L = L 2 ; 

v ; (/xi + ^) 3 r(4 - cw) 

if, on the other hand, 6 = or b = oo, 



mini 



(i o\ J2 '^l i max a min A T 

(6.2) ir = and L = L 

(Hi + ^2)^(4 - a min ) 

where min is the index 1 if b = oo (e.g. if a\ < a 2 ) and is the index 2 if b = 0, max denoting 
the other index. 

3.1 Reflecting Gaussian process as the limit 

Condition 3.1 (lieavy traffic condition,) For each N, let the service rate fi N be given by 

(3.3) n N = NX-f + 9Vn, 
where 6 is some positive constant. 

In addition, we need the below conditions on the distributions of Fx and F 2 : 

(3.4) Fi{x) (2 = 1,2) is absolutely continuous in terms of x; 

(3.5) The density fi(x) (i = 1,2) of Fi satisfies lim fi(x) < 00. 

Before we state our main theorems, we define the scaling processes for each N as follows, 

(3.6) Q N (-) = -j=Q N (-), L N (.)^^=L N (.). 

Theorem 3.1 Under condtions (3.3)-(3.5) and as N — >■ 00, both Q N {■) and L N {■) converge 
in distribution under Skorohod topology to a reflecting Gaussian process Q(-) given by 

(3.7) Q(-) = A( T ) + at (■) - S(X T ) - • +/(■) > 0, 

where the three processes A(j-), S(Xj-) andT(-) are independent each other, and furthermore, 
A(j-) is a Brownian motion with mean zero and variance function Xj-, S(Xj-) is also a 
Brownian motion with mean zero and variance function X^a 2 - , T(-) is a Gaussian process 
with a.s. continuous sample paths, mean zero and stationary increments, whose covariance 
and variance functions satisfy 

(3.8) Cov{f{t),f{s)) = ±(Var(f(t)) + Var(f(s))-Var(f(t-s))), 

(3.9) Var{f{t)) ~ [ ^f* L ® as t ^ 00 for I < a mm < 2, 

[ ir t as t — y 00 and a m i n = 2, 

where H is the Hurst parameter given by H = (3 — a m i n )/2. Moreover, /(•) in (3.7) is a 
non- decreasing process with 1(0) = and satisfies 

POO 

/ Q{s)dl{s) = 0. 
J 
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Remark 3.1 More discussions about reflected Gaussian processes, readers are referred to 
Whitt [29]. From the theorem, we have the following observations. When 1 < a m i n < 2, we 
have that 1/2 < H < 1 which implies that the process T(-) exhibits long range dependence. 
When cti = 2 for i = 1,2, the ON- and OFF-periods both have finite variance and hence 
we have that H = 1/2 and L = 1, which imply that T(-) exhibits short range dependence. 
Finally, the results given in the theorem can be considered as a kind of extension of Little 's 
formula for M/M/l queueing model or considered as satisfying certain state space collapse 
property. 

3.2 Reflecting fractional Brownian motion as the limit 

In this subsection, we suppose that at least one of o\ (i = 1,2) is infinite. To further 
discussion, we need to introduce another time-scaling parameter R and assume that N = 
N(R) goes to infinite as R — > oo. Moreover, we assume that N is taken to satisfy the below 
fast growth condition (and see more discussion in Mikosch et al. [22]) 

(3.10) NRF L (R) -> oo, as R->oo, 

where Fl = F% if L = Li and L is defined in (3.1) and (3.2). Notice that (3.10) implies 

NR 1 ~ olmin L(R) -»• oo. 

Condition 3.2 (heavy traffic condition,) For each N and R, let the service rate /j, R be 
given by 

(3.11) n R = NX-f + 6 (NR l ~ a ™ n L{R)} 1/2 , 
where 9 is some positive constant. 

Next, let da be the normalization sequence given by 

(3.12) d R = (NR 3 - a ™ n L(R)) 1/2 , 
and define 

(3.13) Q R (-) = ^-Q N (R-), L R {.) = ^L N (R-). 

Theorem 3.2 Assuming that conditions (3.10) and (3.11) hold, then as R — > oo, both Q R (-) 
and L R {-) converge in distribution under Skorohod topology to a process Qh{-) given by 

(3.14) Q H (-) = XttB h (-) - 6 ■ +L H (-) > 0, 

where F>h{ ) is a standard FBM, and Ih(-) is a non- decreasing process with 1(0) = and 
satisfies 



roo 

/ Q H {s)dL H (s) = 0. 
J o 
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Remark 3.2 Standard FBM is a mean zero Gaussian process with a.s. continuous sample 
paths and whose covariance structure is as follows 

Cov{B H {t),B H {s)) = \(\t\ 2H + \s\ 2H - \t - s\ 2H ) . 
4 Proof of main theorems 

Let T N (t) be the total cumulative amount of ON time summed over all N sources, that is, 

(4.1) T N {t) ^J2T n (t)= f t W N (s)ds, 

where W N (■) is the superposition of W n (-) for n = 1, N, that is, for each t G [0, oo), 

N 

(4.2) W N (t) = ]T W n (t). 



n=l 



Moreover, let A(t) denote the cumulative number of arrival packets to the server during the 
time interval [0, t], that is, 



(4.3) A{t) = sup { m : ^ «(*) < * f , 



where = 1,2,...} is an exponentially distributed random sequence with mean value 

1/A, which is independent of all processes mentioned before. Then we have the below lemma. 

Lemma 4.1 The stochastic processes A N {-) in (2.5) and A{T N '(•)) in (4.3) have the same 
distribution. 

Proof. To show that A N (■) and A(T N '(•)) have the same distribution, it suffices to show 
that they have the same finite-dimensional distribution for an arbitrary positive integer k 
and arbitrary numbers ti, ■■■,tf ; £ [0, oo) according to Proposition 2.2 in Kallenberg [17]. 

Notice that the process W N (■) in (4.2) takes values in the set N = {0, 1, N} and has 
the piecewise constant sample paths given by 

M 

(4.4) x(t) = < t < G M,n^i / m, 

i=i 

where so, si, sm with so = and sm = oo is a partition of the interval [0, oo) and M is a 
positive integer or infinite. Then we use Dj^[0, oo) to denote the set of all of these functions 
defined in (4.4). Obviously, it is a subset of the Skorohod topological space De[0, oo). Under 
the same topology, Dj\f[0, oo) becomes a measurable space in its own right when endowed 
with the Borel cr-field A n B = {A n B, B G B} where B is the Borel a-field in D E [0, oo) (see, 
for example, Kallenberg [17]). Then there is a probability distribution F N (-) on Dj^[0, oo) 
for the process W N (■) in (4.2), which is uniquely determined by the length distributions of 
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ON- and OFF-periods and the source number N (here, for our purpose, we will not derive 
the explicit expression of F N (-)). 

Basing on the above observation, we first consider the one-dimensional case. For each 
t > and each nonnegative number m, it follows from the independent and stationary 
increment properties of Poisson process that 

P{A N (t) =m}= [ (x> niA (A Si ) = m\W N (-) = x(-)\ F N (dx), 

where x(-) is a sample path as defined in (4.4), N ni \(Asi) is the number of arrival packets 
for the Poisson process with arrival rate n^A during the time interval Asi = min{sj,i} — 
min{sj„i, t} for i G {1, M}, and the integer c is given by c = 1 + sup{i : S{ < t}. Then by 
the independent and stationary increment properties again, we have, 

P{A N (t)=m} = f p{jr NiX (As l )=m\W N (-)=x(-))F N (dx) 

JD M [0,oo) [ i=Q J 

= / piy^NxiiAsi) = m\W N {-) = x(-)\F N (dx) 

Jd m [o,oo) [ i=l J 

= / p{N x (T N (t))=m\W N (-) = x(-)\F N (dx) 

JDj^-10,00) 1 > 

= P{A(T N (t))=m}, 

where Asi is the summation of time intervals during which the arrival rate for the associated 
Poisson process is iX, and r N (t) is the total cumulative amount of ON time from all iV sources 
up to time t along the sample path x(-). 

Secondly, we consider the two-dimensional case (we will omit the discussion for more 
higher-dimensional cases since they are similar). For any t±,t2 G [0, 00) with t\ < t2, and 
nonnegative integers mi and rri2, it follows from the independent and stationary increment 
properties and the definition of conditional probability that 

p{A N (t 1 )=m 1 ,A N (t 2 )=m 2 } 

= f p{A N (t 1 )=m 1 \W N (-) = x(.)\ 

P [A N {t 2 - h) = m 2 - m 1 |VF JV (-) = x(-)} F N (dx) 
= [ P{N x (T(t 1 ))=m 1 \W N (-) = x(-)\ 

P {N x (r(t 2 - ix)) = m 2 - mi \W N {■) = x(-)} F N (dx) 
= p{A(T N (t 1 ))=m 1 ,A(T N (t 2 ))=m 2 }. 

where r{t 2 — t±) is the total cumulative amount of ON time from all N sources during time 
interval [t\,t 2 ) along the path x(-). Hence we have proved that A N (-) and A(T N (-)) have the 
same distribution. □ 
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4.1 Proof of Theorem 3.1 



First of all, we define some scaled and centered processes. For each t > and N > 1, let 



(4.5) 


-N 

A (*) = ■ 


1 

y/N 


(A(JVt) 


- XNt) , 


(4.6) 


S N (t) - - 


1 

s/N 


(s"(t) ■ 


- P N t) , 


(4.7) 


f N {t) ~ ■ 


1 


(r N (t) 


- >yNt) = 



N 



Then we have the following lemma. 

Lemma 4.2 There exist three independent processes A(-), S(Xj-) and T(-) such that 

(4.8) (A N (.),S N (.),f N (.))^(A(.),S(X T ),f(.)) as N oo, 

where A(-) is a Brownian motion with mean and variance function X-, S(X'y-) is a Brownian 
motion with mean zero and variance function X^a^-, T(-) is a Gaussian process with station- 
ary increments, mean 0, stationary increments, whose covariance and variance functions are 
as given in (3.8)-(3.9). 

Proof. First of all, it follows from Functional Central Limit Theorem (e.g., Chen and 
Yao [5]) that 

(4.9) ~A N {.) ee J= (A(N) - XN) => !(•), 



where A(-) is a Brownian motion with mean zero and variance A-. 
Secondly, for each t > 0, we have, 



S N (t) = -L(sup{k:v N (l) + ... + v N (k)<t}- fi N t) 
= -^=(sup{k:v(l) + ... + v(k)<n N t}-fi N t 
= -j=(s?(Nrft)-Ntft), 



'N 

where in the last equation, fi^ is given by 

and Si(-) is the counting process corresponding to the i.i.d. normalized random sequence 
{v(i),i > 1} with mean 1 and variance o~%. It is obvious that ^ — > A7 as N — > 00. Then by 
Functional Central Limit Theorem (e.g., Chen and Yao [5]), we have 

S N (-) S(X r ) as N -»• 00, 
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where S(X'y-) is a Brownian motion with mean zero and variance function X^/cr^-. 

Thirdly, it follows from conditions (3.4)-(3.5) and Corollary 3.1 in Debicki and Pal- 
mowski [11] that the below convergence in distribution is true 



1 N 

(4.10) W N {.) ee -= ]T iW n {.) - jN) 

„=i 

where W"(-) is a stationary centered Gaussian process with a.s. continuous sample paths 
(by Hahn [14] since W n (t) is stochastically continuous) and covariance function ri(-) which 
satisfies (see the proof of Theorem 1 in Taqque et al. [28] for details), 

(4.11) Var Qf* W{u)du^j = 2 j* £ r){u)dudv, 

which has the expression as in (3.9). By Skorohod representation theorem (see, for example, 
Ethier and Kurtz [12]), we can assume that the convergence in (4.10) is u.o.c. Then we have 

f N (-) = ( W N (s)ds -»• / W(s)ds = f (■) u.o.c. as N -»• oo. 

JO JO 

Thus by the definition of weak convergence on C[0,oo) (see, for example, Whitt [29]), Sko- 
rohod representation theorem and Proposition 14.6 in Kallenberg [17], the above u.o.c. con- 
vergence implies weak convergence. Now, we show that T(-) is a Gaussian process. Due to 

(4.11) and Theorem 7 in page 128 of [13], W(-) is mean square integrable in any given finite 
interval [0,T], and therefore it follows from Theorem 3 in page 142 of [13] that T(-) is a 
Gaussian process in [0, T]. Since for any given n € {1, 2, } and any given t\, t n G [0, oo), 
we can find an T\ < oo such that t\, ...,t n belong to the common interval [0, T\\. Hence the 
joint distribution of T(ti), ...,T(t n ) is normal. Thus we can conclude that T(-) is a Gaussian 
process in [0, oo), whose variance function is as shown in (4.11). Since W(-) is stationary, 
T(-) has stationary increments and its covariance function is given by the expression in (3.8) 
due to Proposition 1(b) in Choe and Shroff [6]. 

Finally, by the independence assumptions and definitions of related processes, we know 

-N 

that the three processes A (•), S N (-) and T N (-) are independent each other for each N. Thus 
we can conclude that A(-), S(X'j-) and T(-) are independent each other. Hence we finish the 
proof of the lemma. □ 

To complete the proof of the theorem, for each t > 0, we rewrite (2.8) as the summation 
of centered processes and regulated non-decreasing process as follows, 

(4.12) Q N (t) = X N (t) + I N (t), 
where 

X N (t) = (A N (t) - X 7 Nt) - (S N (B N (t)) - » N B N (t)) - yfNOt, 



I N (t) = fi N f I{Q N (s) = 0}ds. 

JO 
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The process I (•) is non-decreasing process and can increase only when the queue length 
process Q N (■) reaches zero due to the non-idling service discipline and the fact that Q N (t) > 
for all t > 0. 

Lemma 4.3 



X" (•) = -=X» (•) X(-) = A( T ) + AT(-) - S(A 7 -) - (9 • as N -+ oo, 
V 



where A( 7 -) is a Brownian motion with mean and variance function A 7 
Proof. First of all, we prove the following claim to be true 

E N (.) ee J= (^(-) - A 7 7V.) ^ i( 7 .) + AT (■) as N - 

In fact, by Lemma 4.1, it suffices to prove the following claim, 

-j= (A(T N (.)) - A 7 AT.) i( 7 .) + AT (•) as iV 



00. 



00, 



and it is a direct conclusion of Lemma 4.2 and Corollary 13.3.2 of Whitt [29]. Thus, by 
Lemma 4.2 and the independence assumption, we have the below joint weak convergence 



(4.13) 



Moreover, by Skorohod representation theorem, we can assume that the above convergence 
is u.o.c. a.s. Thus it follows from (4.13) that 



(4.14) 



(1^(.), ^(-)) ->(A7-,Ar) u.o.c. a.s. 



Then, due to (4.14), the conditions stated in Theorem 6.5 of [5] are satisfied. So, by the same 
theorem of [5], we know that, for each t > and as TV — > 00, 



(4.15) 



max \B N (s) — s\ 

0<s<t 



0. 



Therefore, by the above discussions and the fact that the associated limiting processes have 
a.s. continuous sample paths, we have 



X N (.)-X(.) 



< 



< 



E 



N 1 



A( T ) + XT(- 



4- s(B N (.))-s(\ T ) 

E N (.)-A( T ) + Xf(. 

S(B N (■))-§ (X T ) 



+ 



+ 



S N (B N (-))-S(XjB N (-)) 
S N (-) - S(X T ) 



+ 



a.s. as N — > 00, 
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where in the second inequality, we used the fact that B N (t) < t for each t > and in the last 
claim, we also used the fact that S(-) is continuous. Thus 

X N (-)->X(-) u.o.c. a.s. as N -»• oo. 

Hence by Proposition 5.3 in Chapter 3 of Ethier and Kurtz [12], the lemma is proved. □ 
Next, similar to the discussion as in (4.12), let 

(4.16) V N (t) = ^=(„ N V N (t)-t), 



and rewrite (2.9) as the summation of centered processes and regulated non-decreasing process 
as follows, 

(4.17) n N L N (t) = Z N {t) + I N (t), 
where 

Z N (t) = ^ N V N (A N (t)) - A N (t)) + (A N (t) - NXjt) - y/NOt. 
Then we have the following lemma. 
Lemma 4.4 

(4.18) (V N (-),A N (-),f N (-)) (V(-),A(X T ),f(.)) as iV ^ oo, 

where V(-), A(-) andT(-) are independent Brownian motions, moreover, V(-) is of mean zero 
and variance function a\-, A(-) andT(-) are given as before. Moreover, 

(4.19) Z N {-) = ^=Z N {-)^X{-) = A{ 1 -) + \f{-)-S{Xr)-e- asN^oo. 

v -/V 



Proof. By applying Functional Central Limit Theorem and the same explanation as in 
Lemma 4.2, one can prove the convergence stated in (4.18). Then it follows from (4.18), 
Lemma 4.1 and random time change theorem that 

Z N {-) = -^Z N {-) V(\ T ) + I( 7 -) + AT(-) -9- as N -»• oo. 
V N 

Notice that V(X'j-) and —S(X^f-) have the same distribution, we can conclude that the claim 
stated in (4.19) is true. □ 

Proof of Theorem 3.1 

Once the above lemmas are obtained, we can go over the following standard procedure to 
finish the proof of the theorem. By Skorohod representation theorem, we suppose that the 
convergence in Lemma 4.3 is u.o.c. Then, by (4.12) and according to Theorem 6.1 in Chen 
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and Yao [5], there uniquely exist a pair of regulated mappings <\> and ip, which are continuous, 
such that for each t > 0, 

I N (t) = -±=I N (t) = J>(x N (t)) = sup (X N (S)Y, 
VA" v 7 0<s<t v 7 

Q N (t) = -^Q N (t) = i>(x N (t))=x N (t) + 4>(x N (t))>o, 

where x~(s) = max{— x(s), 0}. Then by continuous mapping theorem and Lemma 4.3, we 
have, as N — > oo, 

(4.20) I N (-) -»• /(■) = (ir(-)) a.s. u.o.c, 

(4.21) Q N (-) -»• Q(-) = -0 (.*(•)) > a.s. u.o.c. 

Obviously, /(•) and Q(-) have a.s. continuous sample paths, and moreover, /(•) is non- 
decreasing with 1(0) = 0. Since Q N (t) > and I N (t) increases only at times t such that 
Q N (t) = 0, we have for each T > 0, 

(4.22) [ T Q N (t) A ldI N (t) = 0. 



Define 

/ : x G i? ->■ /(s) = x A 1. 

Clearly, we have / G C b (i?). Then by (4.20), (4.21), (4.22) and Lemma 8.3 in Dai and Dai [8], 
we have 

rT 



f Q(t) A ldl(t) = for all T > 0. 
Jo 



Hence /(•) increases only at times t such that Q(t) = 0. 

Finally, by Lemma 4.4 and the same procedure as above, one can prove the weak conver- 
gence for the processes of L N (■) as N — > oo. □ 



4.2 Proof of Theorem 3.2 

Lemma 4.5 Let (3 = 1 — a m i n /2. Then, as T — ^ oo, we have, 

(4.23) U(T) = T^L(TY' 2 -> oo, 

(4.24) V(T) = T ami ™ I 2 ' 1 / 2 / L(Tf/ 2 -»• oo. 

Proof. Since L(T) is a slowly varying function and 1 < a m j n < 2, we know that U(T) is 
a regularly varying function with index < (3 < 1/2, that is, for x > 0, 
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Then, take < e < f3, it follows from Proposition 0.8 in Resnick [25] that there is a fixed To 
such that for x > 1 and T > To, we have 

U(Tx) > {l-e) X P- e U(T). 

Let x — > cc in the above inequality, we know that (4.23) is true. 

Similarly, V(T) is a regularly varying function with index < a m i n /2 — 1/2 < 1/2, then 
by the same reason as above, we know that (4.24) holds. □ 

Now for each t > 0, we rewrite (3.13) as the summation of centered processes and regulated 
non-decreasing process as follows, 

(4.25) Q R {t) = X R (t) + I R (t), 
where 

X R (t) = j- (A N (Rt) - XjNRt) - j- (S N (B N (Rt)) - fi R B N (Rt)) - 9t, 

I R {t) = ^Y N (Rt) = ?£j*I{Q R (s) = 0}ds. 

The process I R ( ) is non-decreasing process and can increase only when the queue length 
process Q R ( ) reaches zero due to the non-idling service discipline. 

Lemma 4.6 For each N and R and under conditions (3.10) and (3.11), we have, as R — > oo ; 
X R (-) converges weakly to a process X(-), that is, 

(4.26) X R (-) =>• X(-) = XttB h (-) - 0- 
where ir and Bh{ ) are given in Theorem 3.2. 

Proof. Due to Lemma 4.1, it suffices to prove the below facts, as R — > oo, 

(4.27) -j-A (t n {R-)) - ^-fi R R ■ 
d R v / dn 

= -i- [A (T N (R-)) - XT N (R-)) + ^-A (t n (R-) - jNR-^j - 9 ■ 

=> \TTB H (-)-6; 

and 

(4.28) ^ (S N (b n (R-)) - n R B N (R-)) => 0. 

As a matter of fact, notice that from the proof of Theorem 1 in Taqqu et al. [28], we know 
that the process T n (-) defined in (2.3) has variance 

Var{T n {t)) ~ 7r 2 t 3 ~ amin L(t) as t -> oo. 
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Then it follows from condition (3.10) and a similar proof as used in justifying Theorem 4 in 
Mikosch et al. [22] that the below weak convergence in the space C[0, oo) is true, 

(4.29) f R (-) = ^-(T N (R-) -jNR-) ^ ttBh(-) as R -»• oo, 

or 

where Bh is standard fractional Brownian motion with H = (3 — a m i n )/2 and 7r is given in 
(3.1) and (3.2). 

Next, by Functional Central Limit Theorem (see, for example, Chen and Yao [5]), we 
have that 

(4.30) A R (.) = j^ji {A(NR-) - XNR-) £«(■) as ^ oo, 

where the weak convergence is in the Skorohod topology and ^ a (-) is a Brownian motion with 
mean zero and variance A. 

Moreover, for each t > 0, we have, 

s R ( t) = ^^(s-m-^Rt) 

= {N l R)1/2 (sup {k : + ... + v(k) < fi R Rt} - fi R Rt) 

= (NRjJl ( S ^ NR ^ ~ NR ^ *) ' 

where in the last equation, /j, R is given by 

/if = A 7 + (R l ~ amm L(R)) 1 / 2 6 , 

and Si (■) is the counting process corresponding to the i.i.d. normalized random sequence 
{f(i),i > 1} with mean 1. Moreover, by Lemma 4.5, we have that /j, r — > A7 as R — > 00. 
Then by Functional Central Limit Theorem, we have 

(4.31) S R (-) £ S (A 7 -) as R -»• 00, 

where £ S (A7-) is a Brownian motion with mean zero and variance Xjcr 2 . 

~~ Ft 

Now notice the independent assumption among the processes A (•), S R (-) and T R (-), 
and the properties that Brownian motion and fractional Brownian motion have a.s. contin- 
uous sample paths, then by Skorohod representation theorem (see, for example, Ethier and 
Kurtz [12]), we can and will assume that the convergence in (4.29)-(4.31) is u.o.c. Thus, by 
Lemma 4.5 and for each t > 0, as R — > 00, 



— > a.s., 
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which implies that as R — > oo, 



(4.32) 



Therefore, by (4.30), (4.32), Random Change of Time Theorem in Billinsley [1] and Lemma 4.5, 
we have 



(4.33) 



^-(A (t n {R-)) - \T N (R-j) 

= U(R)(NR)W (^(^(O)-W(- 



a.s. u.o.c. 



Next, notice that, for each t > 0, 



B?(t) = — ^ < t. 



Then, it follows from (4.31) and Lemma 4.5 that 

j- (S K (B N (Rt)) - » R B N (M)) 



< 



U(R) 
l 

[/"OR) 
a.s. 



(iVfl)V2 



5^ - ji H R- 



Thus, we have, as R — > oo, 

(4.34) j- (S N (b N {R-)) - ii R B N {R-)) -> a.s. u.o.c. 

Hence by (4.29), (4.33) and (4.34), as R -)• oo, the convergence stated in (4.27) and (4.28) is 
true. □ 

The remaining proof of Theorem 3.2 is similar to that used in justifying Theorem 3.1. 
Hence we omit it here. 
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